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Optimizing electronic structure and quantum transport at the 
graphene-Si(lll) interface: An ab-initio density-functional study 
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We use ah initio density functional calculations to determine the interaction of a graphene mono- 
layer with the Si(lll) surface. We found that graphene forms strong bonds to the bare substrate 
and accommodates the 12% lattice mismatch by forming a wavy structure consisting of free-standing 
conductive ridges that are connected by ribbon-shaped regions of graphene, which bond covalently 
to the substrate. We perform quantum transport calculations for different geometries to study 
changes in the transport properties of graphene introduced by the wavy structure and bonding to 
the Si substrate. Our results suggest that wavy graphene combines high mobility along the ridges 
with efficient carrier injection into Si in the contact regions. 

PACS numbers: 73.40.-c, 72.80.Vp, 73.22.Pr, 81.05.ue 



It is now common knowledge that Moore's law, which 
has correctly represented the unprecedented progress of 
Si-based electronics for decades, can no longer be sus- 
tained as device dimensions approach the atomic scale pQ . 
One way to proceed next is to augment Si circuitry 
by taking advantage of the exceptional carrier mobil- 
ity in graphitic nanostructres including graphene or nan- 
otubes [2j [3] . Successful utilization of hybrid devices in- 
volving graphene and silicon necessitates microscopic un- 
derstanding of the morphology, electronic structure and 
transport at the Si-graphene interface. There is reason 
for concern that graphitic carbon may not provide the de- 
sired benefit in this case, since the favorable 7r-bonding 
character has been shown to change in graphene interact- 
ing with related SiC [4HTT] and Si0 2 [12HT4] surfaces. So 
far, only a limited number of studies have investigated 
the interaction between graphene monolayers and pure 
Si. Except for a recent report of successful exfoliation 
of graphene on ultra-clean Si(lll) [15], most studies fo- 
cussed on Si(100) [T6HT9] , where the symmetry difference 
between the overlayer and the substrate raises concerns 
about epitaxy and contact quality. 

Here we study the electronic properties and quantum 
conductance at the graphene-Si(lll) interface. We use ab 
initio density functional theory (DFT) to determine the 
equilibrium morphology of the interface and the nature 
of Si-graphene bonds. We find that the lattice mismatch 
between graphene and Si(lll) can be accommodated by 
buckling the graphene overlayer and creating an array of 
free-standing graphene strips separated by regions cova- 
lently bonded to the substrate. Our ballistic transport 
calculations identify the effect of a covalently connected 
Si substrate on transport in the graphene overlayer and 
describe quantitatively the injection of carriers across the 
interface. 

To gain insight into the equilibrium structure, stabil- 
ity and electronic properties of a graphene monolayer on 
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FIG. 1. (Color online) Optimum geometry and electronic 
structure of wavy graphene on the Si(lll) surface, (a) Equi- 
librium structure of the slab and (b) electron density differ- 
ence An(r) in a plane normal to the surface, (c) Top view 
of the structure. The ridges of C atoms forming parapheny- 
lene chains can be distinguished from sp 3 C atoms covalently 
bonded to Si and C atoms in butadiene like units that are not 
covalently bonded to Si. ai and a2 are the Bravais lattice vec- 
tors defining the 2x1 surface unit cell, (d) Electronic density 
of states (DOS) of wavy graphene/Si(lll) (solid black line), 
wavy graphene only (dashed red line) and planar graphene 
(dotted blue line). E — denotes the position of the Fermi 
level. 



the Si(lll) surface, we performed DFT calculations as 
implemented in the SIESTA code [20] . The surface was 
represented by a periodic array of 6-layer Si(lll) slabs, 
separated by an 8 A thick vacuum region, which were 
connected to a graphene monolayer at the top and ter- 
minated by hydrogen at the bottom, as seen in Fig.Qa). 
We used the Ceperley-Alder [2T exchange-correlation 
functional as parameterized by Perdew and Zunger [22] , 



norm-conserving Troullier- Mart ins pseudopotentials [23] , 
and a double-^ basis including polarization orbitals. The 
reciprocal space was sampled by a fine grid [24] of 
6x12x1 /c-points in the Brillouin zone of the primitive 
surface unit cell and its equivalent for larger supercells. 
We used a mesh cutoff energy of 100 Ry to determine the 
self-consistent charge density, which provided us with a 
precision in total energy of <2 meV/atom. 

Transport properties were investigated using the 
nonequilibrium Green's function (NEGF) approach as 
implemented in the TRAN-SIESTA code [25. Ballistic 
transport calculations for optimized structures were per- 
formed using a single-^ basis with polarization orbitals, 
a 200 Ry mesh cutoff energy, and a 4x60x1 /c-point 
grid [24]. 

Even though silicon and carbon are very similar in 
many ways, graphene is not epitaxial with any silicon 
surface. Previous theoretical studies of graphene on the 
Si(100) surface [18], which has a different symmetry, 
have assumed that the large lattice mismatch may be ac- 
commodated by stretching or compressing laterally the 
graphene overlayer. Since the in-plane compressibility of 
graphene is rather low, the energy cost to enforce epitaxy 
in this way by far exceeds the energy gained by graphene 
bonding to silicon, indicating that graphene should not 
bond to Si(100). 

Also on the Si(lll) surface, which has the same six- 
fold symmetry as the graphene overlayer, there is a large 
11.6% lattice mismatch between the overlayer and the 
substrate. On this substrate, however, there is an alter- 
native way to maintain epitaxy that does not involve in- 
layer compression and still benefits from interface bond- 
ing. When attached to Si(lll), the graphene overlayer 
with the larger lattice constant may buckle and trans- 
form to a super lattice that we call wavy graphene. The 
graphene/Si(lll) superlattice with the smallest 2x1 unit 
cell is shown in Figs. flTa) and flic). We should em- 
phasize that the non-planar, wavy structure of graphene 
in our study is stabilized by strong bonds between sp 3 
hybridized atoms in the overlayer and the substrate, 
which is very different from thermodynamically induced 
rippling observed in graphene on metal substrates [26]. 
Whereas the bare Si(lll) surface is known to undergo a 
7x7 surface reconstruction [27 , no such structural change 
occurs at the graphene-Si(lll) interface, since the dan- 
gling bonds of surface silicon atoms have been saturated 
by forming strong a bonds to the graphene overlayer, as 
seen in Fig. [lib). The rectangular surface unit cell, de- 
limited by the lattice vectors ai and a2, contains 12 C 
atoms in the graphene layer, 12 Si atoms arranged in 6 
slab layers, and 2 terminating H atoms. 

The major benefit of the wavy structure is the coexis- 
tence of ribbon-shaped conducting graphene ridges that 
are detached from the substrate and separated by ribbons 
of carbon atoms bonded to the substrate, enabling car- 
rier injection across the interface. The detached graphene 



ridges contain embedded paraphenylene chains partly re- 
sembling poly-perinaphtalene, and are labeled in this way 
in Fig. [lie). The separating regions contain sp 3 carbon 
atoms covalently connected to the Si substrate and short 
carbon chains resembling butadiene. 

Clearly, changing the period of the wavy graphene 
structure offers a new structural degree of freedom to 
the graphene overlayer. We have investigated the rel- 
ative stability of 2x1, 4x1 and 6x1 supercells of the 
graphene/ Si interface by keeping the bottom four Si lay- 
ers of the slab in the optimum Si bulk geometry [28] . 
Our numerical results allow for a quantitative analy- 
sis of all energy terms associated with the Si-graphene 
bonding. We find that especially the sp 3 carbon atoms 
bind strongly to Si atoms directly underneath, with the 
2.0 A long Si-C bonds comparable to the 1.9 A long co- 
valent bonds in SiC. The covalent bond character is also 
reflected in the electron accumulation in the bond re- 
gion, as seen in the electron density difference An(r) = 
ntot(r) -n grap hene(r) -n S i(in)(r) plotted in Fig.[jjb). If 
we were to attribute the entire graphene-Si interaction 
to these bonds, each of them would contribute 1.62 eV 
towards the binding energy. Obviously, maximizing the 
number of such C-Si bonds is beneficial for the stability 
of the interface. 

To achieve epitaxy, there is an initial energy investment 
associated with the transformation of a free graphene 
monolayer to a wavy graphene structure matching the 
substrate. Even though buckling is less costly than in- 
plane compression, the net energy cost can not be ne- 
glected due the large flexural rigidity and low in-plane 
compressibility of graphene. We find that this energy in- 
vestment decreases with increasing lattice constant ai or 
the corresponding size n of the nxl supercell, favoring 
large supercells [28] . 

The relatively most stable structure of graphene 
bonded to silicon results from an energetic compromise 
between maximizing the number of Si-C bonds and min- 
imizing the buckling energy. Due to the dominant role 
of the strong Si-C bonds, we find that the structure with 
the small 2x1 supercells represents the best energetic 
compromise. 

Graphene will form stable bonds with the silicon sub- 
strate, if the adsorption process is exothermic, i.e. if 

AE = £^tot,graphene/Si — (^tot, graphene + ^tot,Si) < 0- I n 

this expression, ^tot,graphene/Si is the total energy of the 
relaxed wavy graphene structure on Si(lll), £aot, graphene 
that of the equilibrium planar graphene monolayer, and 
-E'totjSi is the total energy of the relaxed Si(lll) surface. 
Defining the average adsorption energy per carbon atom 
as E a d = —AE/Nc, where Nq is the number of carbon 
atoms per unit cell, we find E a d = —0.45 eV in the opti- 
mum case, as the buckling energy dominates over the co- 
valent bonds at the interface. We also found that partial 
hydrogenation of the graphene layer makes the forma- 
tion of a stable graphene superlattice on Si(lll) energet- 



ically much more affordable, as it reduces the adsorption 
energy penalty down to E a & = —0.12 eV in case of 4 
H atoms per C12 unit cell. We expect that additional 
constraints, such as a low density of defects including 
substitutional impurities and vacancies at the interface, 
should turn E a d > 0, yielding a stable bonding geometry 
between graphene and the Si(lll) surface. 

In the following, we will turn to the electronic struc- 
ture and transport in the optimum 2x1 superlattice with 
Ci2Sii2H2 unit cells, shown in Fig. [I] The calculated 
1.56 A corrugation of the wavy graphene normal to the 
surface is sufficient to electronically decouple the car- 
bon atoms in the paraphenylene chains, constituting the 
ridges, from the Si substrate, whereas the remaining car- 
bon atoms in the troughs should be strongly perturbed 
by the vicinity of Si. 

The electronic density of states (DOS) of graphene in 
different environments is shown in Fig. [lid). In compar- 
ison to the free-standing graphene monolayer, which is a 
semi- metal with a smooth DOS near Ep, free-standing 
wavy graphene displays more peaks that reminisce of van 
Hove singularities in ID systems and are caused by a re- 
duction of ppTT interactions normal to the ridges. Apart 
from the only 0.03 eV wide band gap near Ep, the DOS 
of wavy graphene is enhanced with respect to its planar 
counterpart in the ^2 eV wide energy range around Ep 
that is significant for transport. An even larger DOS en- 
hancement near Ep is seen for wavy graphene bonded 
to the Si(lll) surface. Interaction with the substrate in- 
creases the fundamental band gap width to 0.13 eV, in 
analogy to graphene in contact with other semiconductor 
surfaces including SiC [5] and diamond [29]. Results of 
our Mulliken population analysis indicate a small elec- 
tron transfer from silicon to graphene. Such a charge 
redistribution, which is is expected based on the higher 
electronegativity of C as compared to Si, turns the inter- 
face to a pn junction. We find that the extra 0.2 electrons 
per carbon atom are distributed rather evenly across the 
wavy graphene layer. These results all indicate that the 
hybrid graphene/Si(lll) system may display interesting 
quantum transport behavior. 

In order to determine, how contact to a silicon sub- 
strate may affect conduction in a graphene monolayer, 
we performed quantum transport calculations of wavy 
graphene on Si(lll) and present our results in Fig. El We 
distinguished transport normal to the ridges in transport 
geometry A, shown in Figs. |2|a) and^b), from trans- 
port along the ridges in transport geometry B, shown 
in Figs. [2|d) and[2le). We constructed the semi-infinite 
leads of wavy graphene using one cell replicas in geom- 
etry A and two cell replicas in geometry B. The scat- 
tering region consists of three replicas of the 2x1 wavy 
graphene/Si(lll) unit cell, augmented by one additional 
unit cell of wavy graphene on each side to properly de- 
scribe the evanescence of scattering states into the lead 
region. All Si dangling bonds on the surfaces perpendic- 
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FIG. 2. (Color online) Setup for the quantum transport 
calculations for contiguous wavy graphene layers bonded to 
Si(lll). Results for quantum transport normal to the ridges 
in transport geometry A (a-c) are compared to those along 
the ridges in transport geometry B (d-f). (a,d) Schematic ge- 
ometry for the calculations, distinguishing perfect graphene 
leads from the central scattering region, with the direction of 
the current / shown by the arrows. (b,e) Atomic structure of 
the scattering region and its connection to the leads in top 
and side view. (c,f) Quantum conductance G in units of the 
conduction quantum Go as a function of injection energy, with 
E = corresponding to the Fermi level. The conductance is 
given per unit cell normal to the transport direction, shown 
in panels (b) and (e). 



ular to the transport direction have been saturated by H 
atoms. Both leads and the scattering region are infinitely 
wide and periodic normal to the transport direction. 

Transmission spectra G(E) of a contiguous graphene 
layer in different environments are shown in Fig. ^[c) 
for transport geometry A and in Fig. [2lf) for trans- 
port geometry B. In both cases, we compare the quan- 
tum conductance of wavy graphene in contact to Si(lll) 
to that of free-standing wavy or planar graphene mono- 
layers. Our results for geometry A indicate that tr ad- 
mittance normal to the graphene ridges in free-standing 
wavy graphene is reduced to some degree in compari- 
son to planar graphene. The transmission spectrum of 
wavy graphene displays more peaks than that of planar 
graphene, reflecting the changes in the DOS in Fig. [lid) 
including a narrow transport gap of <0.05 eV. Si acts as a 
weak scatterer when connected to wavy graphene. This 
further reduces the conductivity of the wavy graphene 
layer and opens an ^0.35 eV wide transport gap, some- 
what larger than the 0.13 eV wide fundamental band gap 
of the system, seen in Fig. fl[d). 

Electron transmission along the ridges of wavy 
graphene in transport geometry B, shown in Fig. [2Vf ) , 
is greatly enhanced with respect to geometry A. Es- 
pecially impressive is the conductivity enhancement in 
a free-standing wavy graphene monolayer over its free- 
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FIG. 3. (Color online) Geometry and quantum transport cal- 
culations for semi-infinite wavy graphene layers and graphene 
nanoribbons bonded to Si(lll). Results in (a-c) are for trans- 
port geometry C analogous to that of Fig. [5|a), with the 
graphene monolayer disrupted by removing a ridge in the scat- 
tering region. Results in (d-f) are for transport geometry D 
analogous to that of Fig. [5Jd), where removal of every other 
ridge resulted in the formation of bent graphene nanoribbons. 
(a,d) Schematic geometry for the calculations, distinguishing 
free-standing perfect graphene leads from the central scatter- 
ing region, with the direction of the current / shown by the 
arrows. (b,e) Atomic structure of the scattering region and its 
connection to the leads in top and side view. (c,f) Quantum 
conductance G in units of the conduction quantum Go as a 
function of injection energy, with E — corresponding to the 
Fermi level. The conductance is given per unit cell normal to 
the transport direction, shown in panels (b) and (e). 



standing planar counterpart within a broad energy range, 
with the exception of a very narrow transport gap found 
also in geometry A. Even though attachment of the wavy 
graphene monolayer to Si reduces the net conductance of 
the system, this conductance is still higher than that of 
free-standing planar graphene in the ~1 eV wide energy 
window near Ep that is most important for transport. 

Results in Fig. [2] for transport geometry A and B con- 
firm our hypothesis about the formation of anisotropic 
preferential transmission channels in wavy graphene, 
which are responsible for conduction enhancement along 
the conductive ridges containing embedded parapheny- 
lene chains and suppression of conduction normal to these 
ridges. 

To investigate the possibility of charge injection from 
graphene to silicon, we constructed transport geometry 
C by removing a ridge from wavy graphene in the scatter- 
ing region of geometry A, as seen in Figs. |3^a) andJ3^b). 
We followed the approach for geometry A in constructing 
the graphene leads and saturating all Si dangling bonds 
on the surfaces perpendicular to the transport direction 
by H atoms. In absence of the silicon substrate, there is 
obviously no transport due to the gap in a disrupted free- 



standing wavy graphene monolayer. If there were no pos- 
sibility to inject carriers across the silicon-graphene inter- 
face, this would also be true for the disrupted monolayer 
bonded to silicon. Our results in Fig. [3Fc) suggest oth- 
erwise, as we do find transport channels passing through 
the silicon substrate. Obviously, carrier injection across 
the graphene-silicon interface is possible, albeit only into 
and from energetically allowed states below and above 
the 1.1 eV wide fundamental band gap of Si. As in the 
other transport geometries, the transport gap is larger 
than the fundamental band gap. 

Finally, we followed up on our results for geometry 
B, which suggest enhanced conductance along ridges 
of wavy graphene, and studied the effect of laterally 
disconnecting the beneficial paraphenylene-based con- 
ductance channels. Transport geometry D, shown in 
Figs.u3[d) andpfe), has been generated from geometry B 
by removing every other ridge of wavy graphene, creat- 
ing an array of armchair graphene nanoribbons (GNRs) 
that are bent about their axis. These systems have 
been discussed widely as a viable alternative to zero-gap 
graphene [30] 151] . We constructed the graphene leads 
for GNRs and passivated the Si dangling bonds by hy- 
drogen following the approach used for geometry B. Also 
the GNR edges were passivated by hydrogen. 

Our transport results for an array of disconnected arm- 
chair GNRs in transport geometry D are presented in 
Fig. p[f). The reference system, an array of planar 5- 
AGNRs, shows a constant conductance G = IGo corre- 
sponding to one conductance channel in a >2 eV wide en- 
ergy range around Ep, with the exception of an ^0.3 eV 
wide transport gap. These findings agree with previously 
published electronic structure results [30] I5T] . Trans- 
port properties of free-standing 5-AGNRs that are bent 
about their axis are very similar to the planar GNRs in a 
1 — 2 eV wide energy range around the narrow band gap. 
Attaching these GNRs to the Si(lll) substrate causes a 
significant drop in conductance, especially in the conduc- 
tion band region. This result is in stark contrast to the 
related transport geometry B that contains the same con- 
ductive paraphenylene-based ridges as geometry D, but 
does not separate them into nanoribbons. 

The main message of our transport calculation is that 
especially in transport geometry B, the wavy graphene 
monolayer connected to a Si(lll) surface may efficiently 
transport carriers along the graphene-silicon interface. 
Our results for geometry C indicate that the wave- 
function overlap between the overlayer and the sub- 
strate is sufficiently large to permit carrier injection from 
graphene into the valence or conduction band of the sili- 
con substrate. 

In conclusion, we have studied the interaction of a 
graphene monolayer with the Si(lll) surface using ab 
initio density functional calculations. We found that 
graphene forms strong bonds to the bare substrate and 
may accommodate the 12% lattice mismatch by forming 



a wavy structure consisting of free-standing conductive 
ridges that are connected by ribbon-shaped regions of 
graphene, which bond covalently to the substrate. We 
performed quantum transport calculations for different 
geometries to study changes in the transport properties 
of graphene introduced by the wavy structure and bond- 
ing to the Si substrate. Our results suggest that wavy 
graphene combines high mobility along the ridges with 
efficient carrier injection into Si in the contact regions. 
This makes the hybrid graphene-silicon system a suit- 
able candidate for a new generation of high-performance 
electronic circuitry. 
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